      Program Stat
      !*****************************************************************
      Implicit Real*8(A-H,O-Z)
      Parameter(IMax=86)
      Parameter(half=5d-1)
      Dimension rho(IMax),p(IMax),gamma(IMax+1)

      open(3,file='canuto.dat')
      open(1,file='gamma.dat')
      do i=1,IMax
          read(3,*) rho(i),p(i)
      enddo
      do i=2,IMax
          gamma(i)=(log(p(i))-log(p(i-1)))
     +             /(log(rho(i))-log(rho(i-1)))
          write(1,1) half*(rho(i-1)+rho(i)),gamma(i)
      enddo
1     Format(2(1pe12.4))
      close(1)
      stop
      end
